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is  The  correlation  between  the  short  wavelength  features  contained  In  the 
gravity  and  topographic  height  fields  have  been  studied  statistically.  It 
is  found  that  the  slope  of  linear  regression  between  the  two  fields  varies 
with  the  age  of  the  basement  rocks.  This  observation  is  in  qualitative 
agreement  with  the  variations  of  P-wave  velocity  anomalies  with  tectonic 
classes.  These  results  have  been  compared  with  a thermal  model  for  the 
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-^continental  lithosphere.  As  the  lithosphere  moves  away  from  the  ridge 
where  it  is  created  it  thickens  and  induces  the  topography  to  sink  down. 
Assuming  a depth  of  compensation  of  at  most  400  km,  it  is  possible  to 
calculate  the  evolution  of  the  relationship  between  gravity  and  topographic 
height  as  a function  of  age  and  compare  it  with  the  observed  quantities. 

The  geoid  heights  derived  from  the  GEOS-3  altimeter  data  are  also 
studied.  Two  methods  are  being  developed  in  order  to  interpret  the  whole 
range  of  signals  contained  in  the  geoid  heights.  They  each  address  a 
different  class  of  events  and  thus  complement  each  other  in  their  ability 
to  provide  information  on  the  state  of  convection  in  the  earth's  mantle. 

The  long  wavelength  section  of  the  spectrum  yields  information  on  the  depth 
of  the  convection  cells  and  the  viscosity  variations  inside  those  cells 
through  a study  of  the  variations  of  the  admittance  as  a function  of  wave- 
length. The  short  wavelength  section  of  the  spectrum  provides  information 
on  the  time  evolution  of  the  lithosphere,  considered  as  a thin  elastic  plate, 
by  studying  its  response  to  loads  at  several  points  in  its  evolution.  The 
variation  of  the  flexural  rigidity  with  age  will  be  obtained  from  that  study. 
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The  research  performed  under  this  grant  is  to  gain  information  on  the 
lateral  heterogeneities  in  the  earth’s  mantle.  Several  methods  have  been 
employed  for  that  purpose;  they  are  described  in  the  next  section. 


The  first  investigation  was  a study  of  the  correlation  between  the  short- 
wavelength  features  in  the  gravity  field  of  the  earth  and  those  of  the 
topographic  heights.  For  this,  linear  regression  lines  between  free-air 
anomalies  and  topographic  heights  were  obtained  in  5"  x 5”  squares  for  the 
whole  world.  Tne  distribution  of  the  slopes  of  the  regression  lines  was 
then  investigated  statistically  according  to  the  age  of  the  basement  rocks 
for  the  continental  regions  studied.  The  slopes  were  found  to  increase  with 
increasing  age.  That  observation  was  compared  with  models  for  the  continental 
lithosphere  (Crough  and  Thompson,  1976),  in  which  the  heat-flow  equation  was 
solved  to  derive  expressions  for  lithospheric  thickening  and  subsidence  as 
a function  of  age.  Lithospheric  thickening  L'  and  subsidence  S'  can  be 
expressed  as  a function  of  age  t in  dimensionless  form  as  follows: 
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In  the  above,  Q is  the  heat  flux  at  the  base  of  the  lithosphere;  k is  the 
thermal  conductivity;  Tm  is  the  melting  temperature;  p,  p^,  and  pw  are, 
respectively,  the  density  at  temperature  T,  the  density  at  300°K,  and  the 
density  of  water;  Cp  is  the  specific  heat;  and  a is  the  coefficient  of  thermal 
expansion.  We  thus  obtained  the  long-wavelength  behavior  of  the  topography 
and  the  lithospheric  thickness.  Similarly,  using  techniques  discussed  by 
Sclater  and  Francheteau  (1970),  we  obtained  an  expression  for  the  time 
dependence  of  the  temperature,  which  is  written  in  dimensionless  form  as 


V = (1  - z' ) + £ a sin  (mz' ) exp  - (sx')  , (3) 

n 

where 


8„  = (R2  ♦ .V)1'2  - R , 

2 = Lz'  is  the  depth,  x = Lx'  is  the  distance  from  the  ridge  axis  (and  there- 
fore proportional  to  the  age  and  plate  velocity),  T = T T*  is  the  temperature, 
and  R is  the  Reynolds  number.  We  then  derived  the  lateral  density  changes  using 


p = pg(l  - aAT)  . (4) 

Next,  we  evaluated  the  free-air  anomalies  associated  with  columns  that 
are  400  km  deep  (which  is  the  depth  proposed  in  order  for  isostatic  compensa- 
tion to  be  fully  achieved)  and  composed  of  appropriate  proportions  of  crust, 
lithosphere,  and  mantle.  For  a first  approximation,  full  integration  over 
the  whole  plate  seemed  unnecessary.  When  Lambeck  (1972)  performed  such  an 
integration  in  his  study  of  gravity  anomalies  over  ocean  ridges,  he  found 
that  at  a distance  greater  than  200  km  from  a ridge,  the  gravity  anomaly 
associated  with  it  had  decreased  to  insignificance.  Therefore,  we  decided 
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to  calculate  only  those  gravity  anomalies  associated  with  cylindrical  columns 
100  km  in  diameter  and  separated  from  each  other  by  more  than  500  km.  We 
used  the  following  simple  relationship: 

g40Q  = 2itG  x 108  ^PJ400  - l ♦ s'  + (502  + e2)1/2  - 

[502  + (400  + s*)2]1/2f  * U + 50  - (502  + f2)1/2]) 

- 308.6s*  , (5) 

where  6 is  the  gravitational  constant,  pm  and  p^  are  the  mantle  and  litho- 
spheric densities,  respectively  l is  the  lithospheric  thickness  expressed  in 
ki loi>eters,  and  s*  is  the  elevation  in  kilometers;  equation  (5)  works  only 
when  the  continental  regions  that  have  undergone  the  largest  subsidence  are 
taken  as  the  origin.  The  predictions  are  in  qualitative  agreement  with  ob- 
served slopes.  We  further  compared  the  predictions  with  seismic  observations 
using  data  published  by  Sengupta  (1975).  Sengupta’s  work  shows  an  increase 
in  P-wave  velocity  anomalies  with  increasing  age,  in  agreement  with  our 
observations. 

In  the  course  of  the  above  research,  we  became  aware  of  shortcomings  in 
the  data  used,  which  were  1°  x 1°  averages  of  the  gravity  and  topographic- 
height  fields.  Specifically,  the  gravity  data  had  large  uncertainties  and 
poor  coverage  over  large  portions  of  the  oceans,  which  constitute  one  of 
the  easiest  places  to  begin  studying  time  evolution.  Therefore,  we  used 
raw  data,  on  a track-by-track  basis,  from  the  Geos  3 altimeter,  available 
at  SA0.  Utilizing  the  short-wavelength  signals  contained  in  the  geoid 
heights  analyzed  so  far,  we  obtained  information  on  the  time  behavior  of  the 
lithosphere. 

Following  Crough  (1975),  we  can  consider  the  lithosphere  as  a thin  plate 
whose  thickness  increases  with  increasing  time  up  to  a certain  age,  of  the 
order  of  80  m.y.,  and  then  continues  to  increase  at  a progressively  lower 


3 


rate  until  it  reaches  equilibrium  thickness.  Since  the  thickness  of  a plate 
influences  its  mechanical  properties,  it  is  possible  to  study  the  time  evolu- 
tion of  the  lithosphere  by  observing  how  it  deforms  when  loaded  by  seamounts 
placed  at  several  points  along  its  evolutionary  path.  To  examine  the 
mechanical  properties  of  the  lithosphere,  we  assumed  the  thin-plate  model 
developed  by  McKenzie  and  Bowin  (1976).  In  this  model,  the  lithosphere  con- 
sists of  a thin  elastic  plate  overlying  a fluid  medium;  the  plate  is  being 
loaded  by  bathymetric  features  such  as  seamounts,  island  chains,  and  ridges 
and  is  subsequently  deformed.  The  magnitude  and  wavelength  of  the  deformed 
area  depend  mostly  on  the  flexural  rigidity,  which  is  proportional  to  the 
cube  of  the  lithospheric  thickness.  By  studying  the  correlation  function 
between  the  geoid  height  and  the  bathymetry,  we  can  determine  the  flexural 
rigidity  of  the  area  under  investigation.  This  can  be  done  in  one  of  two 
ways:  the  first  is  to  Fourier- transform  the  geoid  height  and  the  bathymetry 
into  wavenumber  space,  divide  the  geoid  height  by  the  bathymetry,  and  obtain 
the  response  function  as  a function  of  wavelength;  the  flexural  rigidity  can 
then  be  deduced  from  the  characteristics  of  this  function.  The  second  method 
is  to  calculate  a theoretical  filter  Z(k)  in  wavenumber  space  by  using  the 
thin-plate  model  (McKenzie  and  Bowin,  1976)  and  varying  the  values  for  the 
flexural  rigidity: 
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In  these  expressions,  pc>  pw»  pm,  and  pg,  are,  respectively,  the  crustal, 
water,  mantle,  and  mean-earth  densities,  r is  the  earth's  equatorial  radius 
t is  the  crustal  thickness,  d is  the  water  depth,  g is  the  average  gravity, 
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F is  the  flexural  rigidity,  n is  the  number  of  points  in  the  filter,  and 
A is  the  spacing  between  consecutive  points  of  the  filter.  The  filter  derived 
in  equation  (6)  is  then  Fourier- transformed  into  direct  space  and  convolved 
with  the  bathymetry,  resulting  in  a theoretical  geoid  height.  The  value  for 
the  flexural  rigidity  that  gives  the  best  agreement  between  predicted  and 
observed  geoid  heights  is  the  one  that  will  be  selected  for  each  area  studied. 

In  practice,  the  method  chosen  will  depend  on  the  type  of  data  available. 
The  first  method  is  more  adequate  when  comparing  gravity  and  bathymetry  data 
from  surface  ships  because  both  sets  of  data  give  equi spaced  points  and  can 
thus  be  easily  Fourier- transformed.  The  second  method,  however,  is  prefer- 
able when  dealing  with  Geos  3 altimeter  data  because  it  is  not  dependent  on 
having  both  bathymetry  and  geoid-height  data  in  a Fourier-transformable  format. 
The  Geos  3 data  are  easily  transformed,  but  the  bathymetry  data  must  be 
reconstructed,  as  rigorously  as  possible,  along  the  subsatellite  position  by 
using  bathymetric  contour  charts;  this  operation  generally  results  in  poor 
accuracy  and  irregular  point  spacing.  The  second  method,  the  two-dimensional 
approach,  is  thus  the  one  we  have  used  to  study  the  evolution  of  the  litho- 
sphere. The  regions  studied  so  far  are  the  Hawaiian-Emperor  Seamount  chain, 
the  Marshall -Gilbert  Island  chain,  and  the  Crozet  Islands,  these  areas  will  be 
part  of  a much  larger  network. 

The  Hawaiian-Emperor  Seamounts  were  selected  first  because  they  have 
been  previously  studied  by  other  methods  from  other  sets  of  data  (Watts  and 
Cochran,  1974;  Walcott,  1976).  This  area  therefore  constitutes  an  ideal 
testing  ground  for  the  two-dimensional  technique.  The  altimeter  passes 
from  Geos  3 selected  in  that  area  are  superposed  on  a map  of  the  region 
in  Figure  1,  and  their  profiles  are  shown  in  Figure  2.  The  profiles  are 
represented  with  respect  to  a reference  geoid  of  degree  and  order  12  calculated 
from  Standard  Earth  IV  spherical-harmonics  coefficients;  they  all  show  the 
features  typical  of  the  region  - a sharp  peak  centered  on  the  island  chain 
flanked  by  a shallow  depression  and  superposed  on  an  asymmetrical  bulge. 
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Figure  1.  Geos  3 passes  studied  in  the  Hawaii an-Emperor  Seamount  region 
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Figure  2.  Observed  geoid-height  profiles  in  the  Hawai i an- Emperor  Seamount 
region,  represented  with  respect  to  a reference  geoid  of  degree 
and  order  12. 
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We  then  calculated  theoretical  filters  using  values  for  flexural 
29  31 

rigidity  ranging  from  10  to  10  dyne-cm;  an  example,  with  a flexural 
30 

rigidity  of  10  dyne-cm,  is  represented  in  Figure  3.  After  convolving 

the  filters  with  the  reconstructed  bathymetry,  we  got  the  results  shown  in 

Figure  4.  The  top  profile  in  the  figure  is  the  observed  geoid  represented 

with  respect  to  a reference  geoid  of  degree  and  order  16,  which  was  chosen 

in  order  to  remove  the  unwanted  long-wavelength  features;  the  middle  profile 

30 

is  the  best-fitting  predicted  geoid,  obtained  with  a value  of  10  dyne-cm 
for  the  flexural  rigidity;  and  the  bottom  profile  is  the  bathymetry  recon- 
structed along  the  subsatellite  positions  from  the  bathymetric  charts 
designed  at  Scripps  Institute  of  Oceanography  by  Chase,  Menard,  and 
Mammerickx  (1970).  Our  values  for  all  passes  are  in  close  agreement  with 
those  determined  by  other  methods  (Watts  and  Cochran,  1974;  Suyenaga,  1977). 

In  the  framework  of  the  time  evolution  of  the  lithosphere,  two  important 
observations  were  made  in  that  region: 

1.  Predictions  for  the  mid-Pacific  Mountains  are  best  when  a lower 
value  for  flexural  rigidity  is  used.  This  is  in  quantitative  agreement 

with  a study  by  Larson  (1.97t>)  on  the  evolution  of  the  Western  Central  Pacific 
Ocean.  Interpreting  the  magnetic  anomalies,  he  characterized  the  mid-Pacific 
Mountains  as  a slowly  spreading  center.  A reduced  plate  thickness  would  thus 
be  expected,  as  observed  in  the  present  work. 

2.  The  flexural  rigidity  associated  with  the  Emperor  Seamounts  is  some- 
what smaller  then  the  easternmost  active  volcanoes  of  the  Hawaiian  chain, 

29 

as  shown  on  Figure  5;  it  is  on  the  order  of  8 * 10  dyne-cm.  The  smaller 
value  can  be  explained  by  taking  into  account  the  age  of  the  seamounts  along 
the  chain  (Clague  and  Jarrard,  1973).  In  a recent  study.  Watts  (1978) 
observed  the  correlation  between  gravity  and  bathymetry  data  obtained  from 
surface  ships  over  several  sections  of  the  Pacific  Ocean:  the  East  Pacific 
Rise,  the  Hawai ian- Emperor  Seamounts,  and  the  Kuril  Rise.  He  deduced  that 
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Theoretical  geoid  filter  calculated  with  a flexural  rigidity  of  10  u dyne-cm. 


Figure  4.  The  top  profile  represents  the  observed  geoid  height  with  respect  to  a reference 

geoid  of  degr°°  and  order  16  in  the  Hawaiian  region;  the  middle  one  is  the  predicted 
geoid,  calculated  with  a filter  of  flexural  rigidity  of  l(pO  dyne-cm;  it  is  con- 
volved with  the  bathymetry  represented  on  the  bottom  profile. 


SCO  «0O  XX)  100  <00  0 >00  100  too  « 0 XX) 


Figure  5.  The  top  profile  represents  the  observed  geoid  with  respect  to  a 
reference  geoid  of  degree  and  order  12  in  the  Emperor  Seamount 
region;  the  middle  profile  is  the  predicted  geoid,  calculated  with 
a filter  of  flexural  rigidity  8 « io29  dyne-cm;  it  is  convolved 
* with  the  bathymetry  represented  on  the  bottom  profile. 
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the  relevant  factor  related  to  flexural  rigidity  is  the  age  of  the  lithosphere 

at  the  time  of  loading.  Although  the  lithosphere  in  the  case  of  the  Emperor 

Seamounts  is  older  than  it  is  at  the  head  of  the  Hawaiian  chain,  the  load 

there  is  proportionally  older,  and  therefore  the  lithosphere  at  the  time  of 

loading  was  younger,  thus  requiring  a smaller  flexural-rigidity  value.  So 

far,  only  one  Geos  3 track  is  available  for  interpretation  in  that  region;  » 

the  difference  observed  is  w'thin  the  range  of  permissible  uncertainties, 

but  we  cannot  draw  any  definitive  conclusions  until  more  tracks  along  the 

seamount  chain  have  been  studied. 

The  two  other  areas  investigated  - the  Marshall  Gilbert  Islands  and  the 

Crozet  Islands  - exhibited  quite  similar  behavior  to  that  found  above.  In 

the  Geos  3 track  shown  crossing  the  Gilbert  Islands  chain  in  Figure  6,  the 

top  profile  represents  the  bathymetry  reconstructed  from  the  Chase  et  al . 

chart,  and  the  bottom  profile  is  the  observed  geoid  height  with  respect  to 

a 16th-order  reference  geoid.  The  Gilbert  Islands  are  the  sharp  feature 

seen  at  1200  km.  Figure  7 shows  a series  of  predicted  geoids  in  the  same 

area  computed  with  several  values  for  flexural  rigidity.  Both  the  width  and 

the  height  of  the  signal  vary  considerably,  and  the  best-fitting  value  in 

30 

this  case  is  between  0.5  and  0.75  x 10  dyne-cm.  The  lithosphere  is  quite 
nld  in  that  area,  of  the  order  of  120  m.y.,  and  therefore  a large  lithospheric 
thickness  is  expected,  which  is  inconsistent  with  the  small  value  found  for 
the  flexural  rigidity.  To  reconcile  the  present  observation  with  Watt's 
model,  it  could  be  speculated  that  the  Gilbert  Islands  constitute  old  loads. 
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The  study  of  the  Crozet  Plateau  was  done  in  collaboration  with  Dr.  Anny 

Cazenave,  from  Groupe  de  Recherche  et  de  Geodesie  Spatiale  and  Centre  National 

d' Etudes  Spatiales,  Toulouse,  who  provided  the  relevant  Geos  3 profiles.  In 

a recent  work,  Cazenave  and  Lambeck  (in  preparation)  analyzed  the  geoid 

anomalies  in  that  region  using  the  three-dimensional  approach  developed  by 

Watts,  Cochran,  and  Selzer  (1975)  for  their  study  of  the  Great  Meteor 

Seamount.  Cazenave  and  Lambeck  found  that  flexural -rigidity  values  ranging 
30 

from  0.75  to  1 x 10  dyne-cm  gave  an  excellent  fit  between  observed  and 
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The  top  profile  represents  the  bathymetry  in  the  Gilbert  Islands  region;  the 
bottom  profile  represents  the  observed  geoid  with  respect  to  a reference  qeoid 

_ r j j » ■»  «-  J 
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Figure  7.  Four  profiles  of  the  predicted  geoid  in  the  Gilbert  Islands  region  calculated 
with  flexural-rigidity  filters  of  0.5,  0.75,  1,  and  2.5  (*  10^0 ) dyne-cm, 
respectively. 
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predicted  geoids.  When  we  applied  the  two-dimensional  approach  described 
above  to  the  Crozet  Islands,  we  obtained  a flexural-rigidity  value  similar 
to  theirs;  this  can  be  seen  by  comparing  the  observed  geoid  plotted  in 
Figure  8 and  the  predicted  geoids  shown  in  Figure  9. 

In  the  Crozet  Plateau  region,  the  age  of  the  load  is  unknown,  and  the 
age  of  the  lithosphere,  according  to  Schlich  (1975),  is  Upper  Cretaceous. 

A comparison  of  Cazenave  and  Lambeck's  method  with  ours  suggests  that  the 
two-dimensional  approach  is  ideal  for  studying  linear  features  such  as  island 
or  seamount  chains,  while  it  offers  less  precision  for  dealing  with  individual 
features.  In  the  case  of  the  Crozet  Plateau,  only  those  tracks  crossing 
the  maximum  altitude  of  the  plateau  gave  a correct  value  for  the  flexural 
rigidity;  all  others  resulted  in  larger  values,  owing  to  the  fact  that  they 
reproduced  only  the  lower  bathymetric  points,  whereas,  in  reality,  the  actual 
observed  qeoid  is  influenced  by  nearby  masses.  In  the  future,  therefore,  the 
feature  being  studied  will  dictate  whether  we  use  the  two-  or  the  three- 
dimensional  approach. 

The  long-wavelength  interpretation  of  geoid  heights  offers  more  complex 
problems.  Several  attempts  were  made  that  did  not  prove  successful.  First, 
the  1°  x 1°  averages  of  the  geoid  height  and  bathymetry  were  compared  directly 
by  using  a program  of  linear  regression.  The  results  of  that  operation  are 
totally  inconclusive,  as  the  comparison  involves  features  belonging  to  a 
whole  spectrum  of  wavelengths,  each  of  which  has  a different  response.  To 
overcome  part  of  that  problem,  we  subtracted  the  long-wavelength  features, 
which  generally  do  not  compare  well  with  bathymetry,  from  the  altimetry  data 
by  removing  a reference  geoid  of  degree  and  order  12.  As  we  had  found  in 
our  work  in  the  Hawaii an-Emperor  Seamount  chain,  this  method  is  effective. 
Comparing  the  residual  geoid  heights  and  the  total  bathymetry,  we  obtained 
short-wavelength  features  in  the  geoid  height  with  a magnitude  of  10  m at 
most;  they  appear  to  be  associated  for  the  most  part  with  distinct  bathymetric 
features. 
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Figure  8.  Observed  geoid  with  respect  to  a reference  geoid  of  degree  and  order  16  in  the 
Crozet  region. 


In  another  effort,  the  1° 


x 1°  averages  for  geoid  heights  and  bathymetry 
were  written  on  Mercator  projections  at  a scale  identical  to  that  on  the 
Mammerickx,  Menard,  and  Sharman  bathymetric  maps  available  from  the 
Massachusetts  Institute  of  Technology  (MIT)  (1.21  per  degree  at  the  equator) 
so  that  we  could  make  a visual  comparison.  This  procedure  also  enabled 
us  to  check  the  reliability  of  the  1°  x 1°  bathymetry  averages  received 
from  Scripps;  they  proved  adequate.  The  geoid  heights  were  then  compared 
with  the  bathymetry  contours  in  the  Hawaiian  region.  We  found  a definite 
correlation  between  shallow  depths  and  higher- than-average  geoid  heights  and 
a general  association  between  large  depths  and  lower- than-average  geoid 
heights.  There  are,  however,  problems  with  that  comparison  - the  relation- 
ship between  geoid  height  and  bathymetry  seems  quite  variable  and  is 
dependent  on  the  background  geoid;  further,  certain  areas  exhibit  very 
marked  local  geoid  features  and  no  sharp  bathymetric  differences.  The  latter 
observation  might  be  caused  by  inadequate  local  coverage,  which  will  be 
rectified  when  the  70  newly  received  altimetry  data  tapes  have  been 
incorporated  in  the  data  ba^e. 

On  the  other  hand,  a more  fruitful  approach  to  the  long-wavelength 
interpretation  of  geoid  heights  would  be  to  study  the  variation  of  the 
admittance  as  a function  of  wavelength,  the  admittance  being  the  ratio  in 
wavenumber  space  of  the  geoid  height  to  bathymetry;  this  technique  has  been 
proposed  by  McKenzie  (1977).  For  that,  purpose,  both  sets  of  data  must  be 
written  in  a coherent  network  and  then  Fourier- transformed.  This  will  be 
the  object  of  future  research. 
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